Sensitivity of the transverse flow towards symmetry energy 
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We study the sensitivity of transverse flow towards symmetry energy in the Fermi energy 
region as well as at high energies. We find that transverse flow is sensitive to symmetry 
energy as well as its density dependence in the Fermi energy region. We also show that the 
transverse flow can address the symmetry energy at densities about twice the saturation 
density, however it shows the insensitivity towards the symmetry energy at densities p/po 
> 2. The mechanism for the sensitivity of transverse flow towards symmetry energy as 
well as its density dependence is also discussed. 
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1 Introduction 



The heavy-ion coUisions (HIC) are the only method presently available in the laboratory to 
produce large volume of excited nuclear matter. The production of such state is essential 
to investigate not only the gross characteristics of nuclear matter but also in exploring 
the explosion mechanism of supernova and cooling rate of neutron stars. After about 
three decades of extensive efforts in both nuclear experiments and theoretical calculations, 
equation of state (EOS) of isospin symmetric matter is well understood by experiments 
of collective flow [1] and subthreshold kaon production [2, 3]. Nowadays, the nuclear EOS 
of asymmetric nuclear matter has attracted a lot of attention. The EOS of asymmetric 
nuclear matter can be described approximately by 

E{p, 5) = E^{p, 5 = 0)+ £;sym(p)5' (1) 

where 5 = is isospin asymmetry, Eo(p, 6) is the energy of pure symmetric nuclear 

matter and Esym(p) is the symmetry energy with Esym(po) — 32 MeV is the symmetry 
energy at normal nuclear matter density. The symmetry energy is E(p,l) - Eo(p,0), ie. the 
difference of the energy per nucleon between pure neutron matter and symmetric nuclear 
matter. The symmetry energy is important not only to the nuclear physics community 
as it sheds light on the structure of radioactive nuclei, reaction dynamics induced by 
rare isotopes but also to astrophysicists since it acts as a probe for understanding the 
evolution of massive stars and the supernova explosion [4]. The existing and upcoming 
radioactive ion beam (RIB) facilities led a way in understanding nuclear symmetry energy. 
Experimentally, symmetry energy is not a directly measurable quantity and has to be 
extracted from observables which are related to symmetry energy. Therefore, a crucial 
task is to find such observables which can shed light on symmetry energy. A large number 
of studies on the symmetry energy of nuclear matter have been done during the past 
decade [5-11]. These studies reveal that in heavy- ion collisions induced by neutron-rich 
nuclei, the effect of nuclear symmetry energy can be studied via the pre equilibrium 
n/p ratio [5-7], isospin fractionation [8,9], n-p differential transverse flow [10,11] and 
so on. These observables have their relative importance depending upon the region of 
density one wants to explore. For e.g. below saturation density (0.3 Po < p < Po), 
observables like fragment yield, isoscaling parameter, isospin diffusion, double n-p ratio 
have been used to extract symmetry energy. On the other extreme, tt'^/tt' ratio, relative 
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and differential collective ffow between triton/He^, n-p differential collective ffow act as 
probe of symmetry energy at high densities. In the low density region, Shetty et al. 
[12] extracted the symmetry energy by comparing the isoscaling parameter from ^°Ar, 
^°Ca+^^Fe, ^^Ni and ^^Fe, ^^Ni+^^Fe, ^^Ni reactions with dynamical Antisymmetrized 
Molecular Dynamics (AMD) calculations [13] and found it to be of the form Egym = 
31. 6(^)''', with 7 = 0.69. Famiano et al. [14] studied the symmetry energy by comparing 
the experimental double neutron to proton ratio in ^^^Sn+^^^Sn and ^^^Sn+^^^Sn reactions 
Isospin-dependent Boltzmann-Uehhng-Uhlenbeck (IBUU) calculations [15] and obtained 
the form Egy^ = 32{-^^)'^, 7 = 0.5. Recently Tsang et al. [7] compared the isospin diffusion 
and n-p double ratio for ^^^Sn+^^^Sn reaction with Isospin-dependent Quantum Molecular 
Dynamics (IQMD) calculations [16] and obtained a similar form of symmetry energy with 
7 = 0.4-1.05. The situation is worse at higher densities. The results are model dependent 
and contradicting also. The FOPI collaboration at GSI studied tt^/tt" ratio in ^°Ca-|-^°Ca, 
'^^Ru+'^^Ru, '^^Zr+^eZr, and ^^^Au+^^^^Au reactions [17]. A comparison of this data [18] 
with IBUU calculations showed a softer form of the density dependence of symmetry 
energy, which is in contrast to those obtained from the low density studies where stiffer 
form of symmetry energy reproduced the data well, tt'^ /Tr~ ratio was also compared with 
Improved Quantum Molecular Dynamics (ImQMD) calculations by Feng et al. [19] and 
favored a stiffer form of symmetry energy and in contradiction with IBUU calculations. 
At densities higher than saturation density, collective flow has also been proposed as a 
novel mean to probe the high density behavior of symmetry energy [10]. In this paper, we 
aim to see the sensitivity of collective transverse in-plane flow towards symmetry energy 
and also to see the effect of different density dependence of symmetry energy on the 
same. The various forms of symmetry energy used in present study are: Egym oc Fi{u), 
Egym oc F2{u), and Egym oc F3{u), where u — Fi{u) oc u, F2{u) oc u'^'^, F^lu) oc u^, and 
F4 represents calculations without symmetry energy. The different density dependences 
of symmetry energy are shown in fig. 1. The various lines are explained in the caption 
of the fig. 1. Section 2 describes the model in brief. Section 3 explains the results and 
discussion and Sec. 4 summarizes the results. 
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Figure 1: (Color online) The density dependence of symmetry energy for various forms of 
symmetry energy : Fi (7 = 1.0) (solid), F2 (7 = 0.4) (dashed), and F3 (7 = 2.0) (dotted). 



2 The model 

The present study is carried out within the framework of isospin-dependent quantum 
molecular dynamics (IQMD) model [20]. The IQMD model treats different charge states 
of nucleons, deltas, and pions explicitly, as inherited from the Vlasov-Uehling-Uhlenbeck 
(VUU) model. The IQMD model has been used successfully for the analysis of a large 
number of observables from low to relativistic energies. The isospin degree of freedom 
enters into the calculations via symmetry potential, cross sections, and Coulomb interac- 
tion. 

In this model, baryons are represented by Gaussian-shaped density distributions 

1 1 2L 

fii^'^P^t) = ^^^P^~[^"'"'^(^)]'^) ^ exp(-[p-pl(t)]2— ) (2) 

Nucleons are initialized in a sphere with radius R = 1.12 A^/^ fm, in accordance with 
liquid-drop model. Each nucleon occupies a volume of h^, so that phase space is uniformly 
filled. The initial momenta are randomly chosen between and Fermi momentum (pf)- 
The nucleons of the target and projectile interact by two- and three-body Skyrme forces, 
Yukawa potential, and Coulomb interactions. In addition to the use of explicit charge 
states of all baryons and mesons, a symmetry potential between protons and neutrons 
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corresponding to the Bethe-Weizsacker mass formula has been included. The hadrons 
propagate using Hamilton equations of motion: 

dfi d{H) dpi d{H) 



dt dpi ' dt dfi 

with 



(3) 



(H) = {T) + {V) 

- E^ + EE / Mr,P,t)VHf',f) 

i i j>i 

X fj {r',p',t) df df' dpdp'. (4) 

The baryon potential V^j, in the above relation, reads as 

V ^{f — r) = Vg^yj.jjjg + Vy^]^^^ + V'j^o^j + Vgy^ 

= M(f'-r-)+M(r'-rV-^(^^)] 



+^3 



exp(|(r ' - r)|///) Z,Z/: 



,2 



(|(f'-f)|//.) \{f'-f)\ 



+U-T3iT3jS{ri' -fj). (5) 

Here Zi and Zj denote the charges of ith and jth baryon, and T^i and T^j are their 
respective components (i.e., 1/2 for protons and —1/2 for neutrons). The parameters fi 
and ti,...,t4 are adjusted to the real part of the nucleonic optical potential. For the density 
dependence of the nucleon optical potential, standard Skyrme type parametrization is 
employed. 

3 Results and discussion 

We simulate several thousands events for the neutron-rich system of ^^Ca+^^Ca and 
™Ca-|-^°Ca at energies of 100, 400, and 800 MeV/nucleon at impact parameter of 
b/bjnax=0.2-0.4. Wc use a soft equation of state and isospin- and energy-dependent cross 
section reduced by 20%, i.e. a = 0.8 cr^nn^. The details about the clastic and inelastic 
cross sections for proton-proton and proton- neutron collisions can be found in [20, 21]. The 
cross section for neutron-neutron collisions is assumed to be equal to the proton-proton 
cross section. 

Since ^°Ca has a very high asymmetry, so to ensure the stability of the nuclei in the 
present study, we display in fig. 2 the time evolution of root mean square radius of single 
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Figure 2: (Color online) Time evolution of root mean square radius of a single ^°Ca, ^^Ca 
and ^°Ca nuclei in coordinate (top panel) and momentum space (bottom) obtained with 
IQMD for EOS used in the present study for symmetry energy Fi(u). 



nucleus of ^°Ca (sohd line), ''^Ca (dashed), and ^°Ca (dotted) in the coordinate (fig. 2(a)) 
and momentum space (fig. 2(b)). The results are displayed for nuclei intialized with 
symmetry energy Fi(u). We find that the stability is of same order for all the three 
nuclei. 

There are several methods used in the literature to define the nuclear transverse in- 
plane flow. In most of the studies, one uses (p^M) plots where one plots (px/A) as 
a function of Vc.m./^beam- Usiug a linear fit to the slope, one can define the so-called 
reduced flow (F). Alternatively, one can also use a more integrated quantity "directed 
transverse in-plane flow {p'^")" which is defined as [20]: 

= ^E^ig^^^(^)>^-(*)' (6) 

i 

where Y{i) and Px{i) are the rapidity distribution and transverse momentum of the ith 
particle. In this definition, all rapidity bins are taken into account. It, therefore, presents 
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Figure 3: (Color online) Left panel: < ^ > as a function of Yc.m./^beam for different 
energies of 100, 400, and 800 MeV/nucleon for different forms of symmetry energy. Right 
panel: The time evolution of < jf^^ > at 100, 400, and 800 MeV/nucleon for different 
forms of symmetry energy at b/bmax=0.2-0.4. Lines are explained in the text. 



an easier way of measuring the in-plane flow rather than complicated functions such as 
(Px/A) plots. 

In fig. 3 we display the < ^ > as a function of Yc.m./^beam at final time (left panels) 
and time evolution of < p'^^^ > (right panels) calculated at 100 (top panel), 400 (middle) 
and 800 MeV/nucleon (bottom) for different density dependence of symmetry energy. 
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Figure 4: (Color online) The time evolution of rapidity distribution for the calculations 
with no symmetry energy for various bins at b/bmax=0.2-0.4. Lines are explained in the 
text. 



Solid, dash dotted, and dotted lines represent the symmetry energy proportional to p, p°'^ 
and whereas dashed line represents calculations without symmetry energy. Comparing 
left and right panels in fig. 3, we find that both the methods show similar behavior to 
symmetry energy. For eg. at incident energy of 100 MeV/nucleon for Esym oc p^ "^, < p'^" > 
= 0. Similarly, the slope of < ^ > at midrapidity is zero. We also find that the transverse 
momentum is sensitive to symmetry energy as well as its density dependence Fi(m), F2{u) 
and F^lu) in the low energy region (100 MeV/nucleon). At energies above Fermi energy, 
both the methods show insensitivity to the different sjTiimetry energies. This is because 
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the repulsive nn scattering dominates the mean field at high energies. 

To understand the sensitivity of transverse momentum to the symmetry energy as well 
as its density dependence in the Fermi energy region, we calculate the transverse flow as 
well as rapidity distribution of particles having < 1 (denoted as bin 1) and particles 
having > 1 (bin 2) separately at all the time steps. Since both the methods show 
similar behavior to symmetry energy, so for simplicity the following discussions will be in 
terms of < p^" >. 

In flg. 4 we display the rapidity distributions at 100 MeV/nucleon of all the particles 
(dotted line), particles corresponding to bin 1 (solid) and to bin 2 (dashed) at 0, 10, 20, 30, 
40, and 60 fm/c. We have calculated rapidity distributions for different forms of symmetry 
energy used in this paper. We find that it is insensitive towards the symmetry energy 
[22, 23]. In fig. 4, we display the rapidity distribution calculated without symmetry energy. 
During the initial stages we see the two Gaussians at projectile and target rapidities for 
all the 3 bins. The peaks of Gaussians will be more prominent at higher energies. The 
interest for our discussion is in bin 1 and bin 2. During the start of the reaction (0 fm/c) 
more number of particles lie in bin 1 i.e. more number of particles have ^ < 1. As the 
nuclei begin to overlap, the density increases in the overlap zone. Now, the number of 
particles increases in bin 2 (at 10 fm/c). Prom 10 fm/c to 20 fm/c, the number of particles 
keep on increasing in bin 2 at midrapidity, i.e particles from large rapidity keep on shifting 
to the bin 2 in the midrapidity region. This is expected since at incident energies in the 
Fermi energy region dynamics is governed by the attractive mean field. The dominance 
of attractive mean field will prompt the deflection of particles into negative angles i.e. 
towards participant zone. After 20-30 fm/c, the expansion phase of the reaction begins 
and number of particles keep on increasing in bin 1 and by 60 fm/c most of the particles 
lie in bin 1. 

In flg. 5 we display the time evolution of < p'^^^ > for different symmetry energies 
used in this paper at 100 MeV/nucleon for particles lying in the different bins. Lines 
have the same meaning as in fig. 4. Panels a, b, and c are for Egym oc p,p^'^, and p^, 
respectively. Panel d is for calculations without symmetry energy. The total < > 
is negative during the initial stages and keep on decreasing till 30 fm/c which indicates 
dominance of attractive interaction. In panels a and b, it becomes positive whereas in 
panel c and d it remains negative during the course of the reaction. If we look at < p^"^ > 
of particles lying in bin 1 for Fi{u) (flg. 5a) and F2{u) (flg. 5b) in the time interval to 



9 




20 40 60 80 20 40 60 80 
Time (fm/c) 

Figure 5: (Color online) The time evolution of < p'^" > for different forms of symmetry 
energy for different bins at b/bmax=0.2-0.4 . Lines have same meaning as in fig. 4. 

about 20-25 fm/c, we see that it remains positive. It increases with time upto 15 fm/c 
and reaches a peak value. This is because in the spectator region (where high rapidity 
particles lies) the repulsive symmetry energy will accelerate the particles away from the 
overlap zone in the transverse direction. After 15 fm/c, < p^"' > (of particles in bin 1) 
begins to decrease. This is because these particles will now be attracted towards the 
central dense zone. As shown in fig. 4, from 10 to 20 fm/c particles in bin 2 keep on 
increasing in the midrapidity region. In case of Fi{u) and F2(m), particles which enter the 
central dense zone (bin 2) have already a high positive value of < p^''' > (i.e. going away 
from dense zone). So attractive mean field have to decelerate the particles first, make 
them stop and then accelerate the particles back towards the overlap zone. At about 
20-25 fm/c particles from bin 1 have zero < p'^^^ > (see shaded area in fig. 5a and 5b). Up 
to 30 fm/c, particles feel the attractive mean field potential after which the high density 
phase is over, i.e. in case of Fi(u) and F2{u) between 0-30 fm/c particles from bin 1 
are accelerated towards the overlap zone only for a short time interval of about 5 fm/c, 
whereas for the case of F^lu) (fig. 5c) and F4 (fig. 5d) between 0-30 fm/c, particles from 
bin 1 are accelerated towards the overlap zone for a longer time interval of about 20 fm/c 
between 10-30 fm/c. Moreover, the < p'^^^ > of particles lying in the bin 1 (for F3(m) and 
F4) follows a similar trend. This is because, for p/po < 1 the strength of symmetry energy 
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Figure 6: (Color online) (a) The time evolution of < p'^^^ > at 100 MeV/nucleon for 
different forms of symmetry energy for ^'^Ca+'^^Ca. (b) The time evolution of < Pmax/ Po > 
and (c) < Pavg/po > for Fi(''^) and F4. Lines have same meaning as in fig. 3. 

F3(-u) will be small and so there will be less effect of symmetry energy on the particles 
which is evident from fig. 5c where one sees that the < p'^" > remains about zero during 
the initial stages between zero to about 10 fm/c. 

The < p^'^ > due to particles in bin 2 (dashed line) decreases in a very similar manner 
for all the four different symmetry energies between 0-10 fm/c. Between 10-25 fm/c, 
< p'^" > for F3(m) and F4 decreases more sharply as compared to in case of Fi(u) and 
F2(m). This is because in this time interval particles from bin 1 enters into bin 2. As 
discussed earlier, < p'^^^ > of particles entering bin 2 from bin 1 in case of Fi{u) and F2(m) 
will be less negative due to stronger effect of symmetry energy as compared to in case of 
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F^iu) and F4. 

Since the reaction ®°Ca+®°Ca is an extreme case with large isospin asymmetry. To 
check if the above predicted effects survive in reactions which are experimentally accessi- 
ble, we simulate the reaction of ^^Ca+^^Ca for all the different forms of symmetry energy 
used in the present study. The reaction ^^Ca+^^Ca has been used in many previous studies 
[24]. The results are shown in fig. 6(a). We find that even for this reaction the transverse 
flow shows sensitivity to symmetry energy as well as its density dependence. In fig. 6(b) 
and 6(c), we display the time evolution of maximum density < Pmax/po > and average 
density < Pavg/po >, respectively at 100 MeV/nucleon. We find that the density reached 
is about 2.0 times the saturation density. Moreover, the maximal density is reached in 
the time interval 0-30 fm/c and the effect of symmetry energy on < p'^" > of particles 
during this interval decides the fate of the final value of < p'^^^ >. Thus transverse flow 
can address the symmetry energy at densities about 2.0 times than saturation density. 

4 Summary 

We have checked the sensitivity of transverse flow towards symmetry energy in the Fermi 
energy as well as at high energies. We have found that transverse flow is sensitive to 
symmetry energy as well as its density dependence in the Fermi energy region. We have 
also shown that the transverse flow can address the symmetry energy at densities about 
twice the saturation density, however it shows the insensitivity towards the symmetry 
energy at densities p/po > 2. We have also discussed the mechanism for the sensitivity of 
transverse fiow towards symmetry energy as well as its density dependence. 

This work has been supported by a grant from Centre of Scientific and Industrial 
Research (CSIR), Govt, of India and Indo-Frcnch Centre For The Promotion Of Advanced 
Research (IFCPAR) under project no. 4104-1. 
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